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1 Introduction 



Recently, the universality of the wide distance two level connected correlator in random matrix the- 
ories is shown by Ambj0rn, Jurkiewicz and Makeenko Brezin and Zee regarded this universality 
as importance from the viewpoint of physics in disordered systems . Since the correlator does not 
depend on the probability distribution of a matrix ensemble, we may use the correlator calculated 
in the corresponding simple Gaussian ensemble with the same symmetry for the non-trivial disor- 
dered system in concern. Compared to the short distance correlator which is already well-known as 
a universal quantity 0], the wide correlators are able to be calculated explicitly in extensive types 
of ensembles in various ways ^. Actually, one observes the strongly universal properties which 
random matrix theories have themselves in the wide connected correlators. Their universality classi- 
fication is also done in 0. These mathematical studies of random matrix theories themselves enable 
us to recognize the real universal nature of level statistics. 

In this paper, we examine the replica method to calculate the wide connected two level correlator 
in gaussian orthogonal (GOE), unitary (GUE) and symplectic (GSE) ensembles. These simple en- 
sembles GOE, GUE and GSE can describe a time reversal system without spin, a general system and 
a time reversal system with spin, respectively |^ . The replica method is well-known as a convenient 
scheme to calculate two body Green's function in some models of the Anderson localization Even 
though this method cannot work well to calculate the short distance correlator ^j, here we show this 
method is useful to calculate the wide two level correlator in the expansion. This method is 
much simpler than the supersymmetry method which enables us to calculate both the wide and short 
distance correlators Here, we remark the reason why the replica method is not good for the 

short distance correlator. We are going to calculate the explicit form of the two-level and three-level 
correlators in GOE, GUE and GSE, which are identical to those calculated by solving functional 



equations |TT| or by diagrammatic method ||T2|. Here we study a matrix theory defined by the 
following probability distribution 

P{H) = ^ exp (-f Tri/^^ , Zh ^ J DH exp (-f , (1) 

where H is a.n N xN matrix. The explicit form of the measure DH depends on a type of an ensemble 
of H. Our interest is focused on computing averaged Green's functions by the replica method. We 
show their calculation method explicitly for GOE in section 2, for GUE in section 3 and GSE in 
section 4. 
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2 Gaussian orthogonal ensemble 

To begin with, we explore the case of the GOE which is defined by the ensemble of real symmetric 
matrices obeying the probability distribution ([l|). The measure DH is explicitly written as 

N 

DH =X{dH^kX{dHi,. (2) 

k=l i<j 

2.1 One-point function 

We calculate the averaged one-point function 

Note that —ie{e> 0) means that G{z) considered here is the advanced Green's function. We assume 

1^1 > V2 (4) 

in the case of the GOE. The meaning of this assumption will be clarified later. Results for \z\ < ^J2 
will be extracted employing the analytic continuation that is uniquely determined by —ie. 

In order to apply the replica method, we first introduce the n-fiavor real vectors having N 
component (0")i<i<Ar, where the superscript a, which runs from 1 to n, specifies the fiavor. Let 
us consider the following integration over for constructing a generating function of Gi^zy. 

nz,H]^^^l D<P exp{-^ *r {z-H- te) ^ }, Z^ = j Dct> exp (-1 , (5) 

where * means transposition and the summation over the repeated indices are implicitly taken. The 
measure is defined by 

n N 

= n n ^c- (6) 

a=l t=l 

Note that the integration over for F converges owing to —it. Hereafter we do not explicitly write 
—ie for brevity. 

Taking the derivative with respect to 2;, we have 

^^1^ = ^det-/^(,:(. - ff)) 

OZ OZ 

= -^Tr-^det-"/2(.(^ -//)). (7) 
We define W{z) as the average of F over H: 

W{z)^{T[z,H]). (8) 



2 



It follows from eq. (j^) that 

«(^)=i";j(-j7)i'^"w- w 

which indicates that W{z) plays a role of a generating function of G{z). Note that we need to take 
the limit n — in order to eliminate the determinant factor in eq. (^. 

Let us compute W{z). The average over H in eq. (§) is explicitly carried out and we have 

w{z) = yJ^^ {-^Yzr - ^( W)') • (10) 

The quartic interaction can be represented as 

exp (-^( W)') = ^ / DQexp L^Q,,Q,, + '_Y^bQ\ ^ ^^^^ 

where {Q)ab is the real symmetric n x n matrix and the normalization factor is defined by 

Zq = J DQexp l^-^QabQb?j . (12) 



1 /■ / N i \ 

W{z) = J DQD<P exp (-yTrQ^ - - g)0,j . (13) 



Using the representation (0) we get 

1 /■ „ / N, 

Note that indices a and b are implicitly summed in the above notation: we regard (pi as an n- 
component vector. 

The 0-integral becomes Gaussian and we find 



Wiz) = ^JdQ exp{-^Tr(logz(;. - Q) + Q')}. (14) 

Q 

It can be written as integration over eigenvalues |^ 

^(^) = ^lf{d^a \Hu)\ exp{-iVX:^7(tXa)}, (15) 

a=l a=l 

where C„ is a constant that depends on n and A(m) is the Van der Monde determinant 

A(m) = \{{Ua-U,). (16) 
a<b 

The function g appeared in the exponent is defined by 

gi.^) = ]j^{'^ogi{z - x) + x^). (17) 
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Now we compute W{z) by the 1/A^-expansion. Let us first compute it in the leading order. Since 



the contribution from the Van der Monde determinant in eq. (^) is neglected in the leading order 



Ma's are completely decoupled each other and W{z) is evaluated by a single- variable integration, i.e., 

W{z) ~ const. (^J du exp {-Ng{u)}^ (18) 

in the leading order. We can evaluate the right-hand side of ([18|) by the saddle-point method. The 
saddle-point equation is 



We have the two solutions m±, where 
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u± = -{z±^z^-2). (20) 



We find that contribution to the integral ([18| ) from u ^ dominates over that from u ~ In 
fact, a straightforward calculation gives 

Re g{u-) < Re g{u+), if z > V2. (21) 

Therefore we get 

W{z) ~ const. e-^"3("-) (22) 
in the leading order. Inserting the above result into eq. (^, we obtain the well-known result 



G{z) = z- vG^^, (23) 

where the overall constant in eq. ( p2D is chosen such that G{z) should satisfy the boundary condition 
G{z) ^ 1/z as 2; — > cxD. 

We can compute higher-order corrections by expanding g{ua) in eq. (K^ around u = U-. 



1 °° AI~"/2 

g{ua) = g{u.) + ^g"{u.) yl + Y. ^9^''\u.) y^, (24) 

n=3 ^• 

where we denote 

Ua-u^=N-^'^ya. (25) 

Here we shall calculate up to the first-order correction, which is sufficient for determining the 
connected two-point function, as we will see below. To this end, we can neglect the cubic and the 
higher terms in eq. (^4]). Then we find 



r 1 " 

ly(^) ~ C„e-^"^("-) / \[dya \^{y)\eM--Y.9'\u-)yl} (26) 

•' " ^ a=l 
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up to the next-to-leading order. Note that 

= ^fW^(z-^/7^^2), (27) 

which is real and positive under the assumption (^). Therefore the integration over is easily 
performed going back to representation of an integration over real symmetric matrices. Since there 
are n(n + l)/2 independent variables in a real symmetric matrix, we get 

" yn^2/JA(?/)|e^P{-2E/(^-)^a} = ^y^?exp{--/(n_)^g„M 



'^Q ■> a ^ a=l •> ^ a,b 

/ ^ \n(n+l)/4 



:28) 



Thus we conclude 

/ , X n(n+l)/4 

Wiz) ~ e-^-("-) (^^j (29) 

up to the next-to-leading order. Using eqs. (^ and (P?]), we obtain G{z) including the first-order 
correction: 

Giz) = {z- v/^^) (l + + 0{1/N'). (30) 



This result is identical to that obtained in 11 



Here we comment on the validity of the saddle point method in eq. (^). The saddle point 
evaluation in the 1/A^ expansion holds only in the region where g"{uJ) is 0{N^). This approximation 
becomes incorrect for z = ±\/2 + 0{1/N) which yields g"{uJ) = 0{\/N). 



2.2 Two-point function 

Now we turn to the two-point function 

We shall show how to compute wide correlation by the replica method, i.e. , we assume zi — Z2 ^ 
0{N^) as well as the condition @j. 

As is the case of the one-point function, we start with the following integral 

r[zi, Z2, H] ^ -i— / D0iD02 exp{-^ Vi (^1 -H)ri-'i: %l (^2 - H) 0^}. (32) 
^4>i^4>2 J ^ ^ 

Here we have introduced the two species of vectors labeled by 1 and 2. The two species respectively 

have n- and m- flavor, which means that 1 < a < n and that 1 < p < m. We average Z[zi^ z^^ li\ 

over H and deflne W{zi,Z2), a generating function of G{zi, Z'l). 

W{z^,Z2) = i^{zx,Z2.E\). (33) 
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The two-point function is derived from W{zi, Z2) in the following formula: 



n,m~ 



Our strategy for computing W{zi, Z2) is the same as in the case of the one-point function. Namely, 
we first perform the integration over H of W{zi,Z2). Second, we introduce the auxiliary matrix Q 
in order to carry out the integration over 0i and 02- The integration over the vector variables bring 
us the theory described by Q, which can be investigated by the saddle-point method for large N. 

The first step is easily performed. We have 

W{z,,Z2) = ^^jD<P^D<P2exp{-'- (35) 
Next, let us introduce the (n + m) x (n + m) real symmetric matrix Q using the following expression: 



e H j: i . (a^) 



where Q^^ and Q^^ are nxn and mxm real symmetric matrix respectively. Since Q is real-symmetric, 
Q^^, which is an n X m matrix, satisfies Q^^ = In order to write down a formula corresponding 
to (0), we employ the following notations: 



- I 1 , (37) 



which is (n + m)-component vectors, and 

° I, (3^ 

Z2l„ ' 



where /„ is the unit nxn matrix. We can check that the right-hand side of eq. (pSD is identical with 

1 . iV. 



J DQD4>,Dc^2 exp{--Tr - (2 - g) (39) 

After integrating the vector variables, we acquire 

W{z^, ^2) = ^ / DQ exp |-^Tr (log*(2 - g) + Q^^j . (40) 

We analyze the above integral formula by the saddle-point method. The saddle-point equation is 

{z-Qy'=2Q. (41) 
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Or equivalently, 

2g"(^i-g")-2gi2g2i ^ ^ 

-2Q^^Q^^ + 2Q^\z,-Q'^) = 

-2Q^'Q'' + 2Q'\z2-Q^') = 1. (42) 

After transposing the third equation and subtracting it from the second equation, we find 

Qi2 ^ q21 ^ Q_ (43) 

The remaining equations containing Q^^ and Q^^ are solved by diagonahzation. Suppose that Q^^ 
and Q"^^ are respectively diagonalized by Oi and O2. Explicitly, 

*0i Oi = diag (Ml, *02 g'' O2 = diag (wi, ■ ■ ■ , t;^). (44) 

Then the equations to be solved are 

1 



2{Zi - Ua)' 

2(2:2 - Vp) 

which are the same form as the saddle-point equation appeared in the previous subsection. Following 
the argument on the most dominant saddle-point that is carried out in the one-point function, we 
should choose the saddle point 

Ua = H-, for all a, 

Vp = V-, for all p. (46) 



Performing the inverse transformation of eq. (^4]) , we obtain the most dominant saddle point Q in 
the original basis: 

Q=f"-/" ° ]. (47) 

Next we consider fluctuations around the saddle point Q. Using the identity 

Tr log(A + 6 A) = Tr log A + TtA-^6A - ^TtA-^6A A'^6A + ■■■, (48) 
and the saddle point equation (^), we have 

Tr (^\ogi{z - Q - i^) + (Q + i^)2 j = Tr (log z(2 - Q) + Q') + ^Tr -2SQ Q 6Q Q) + ■ ■ ■ . 

(49) 



7 



Inserting the explicit form ( |T7| ) into the above, we see that the generating function becomes 
W{zi,Z2) ~ ±e-^"^("-) / D{5Q^^) exp{-^(l - 2uJ)Ti f^g")'} x 
e-''"'3{v^^ I D{5Q^^) exp{-i(l -2t;_2)Tr ((5g22)'} x 
j D{6Q^^) exp{-(l - 2u_v.)Tt {SQ^^f} (50) 

up to the next-to-leading order. Remembering that g"{uJ) = 1 — , we get 

\ n(n+l)/4 / \ m(m+l)/4 / i \ "ni/2 

\ nm/2 

^_2u v ) • ^^^^ 

Inserting the above formula into eq. (0), we can derive the two-point function. Differentiating W{zi) 
and W{z2) brings about the disconnected part, G{zi)G{z2), while the last factor contributes to the 
connected part. Thus we obtain 

G{z,, Z2) = G{z,)Giz2) - Jj^^ log (l - lG{z^)Giz2)) , (52) 
where we have used eq. (||) and 2m_ = G{zi) + 0{1/N). This result with respect to the connected 



part agrees with that obtained by solving functional equations H, 11 1 and by a diagrammatic method 



12| . Especially the disconnected part up to order in eq. (|5^ is consistent with ||TT[ . 



The saddle point evaluation holds only in the case of 1 — 2'u_f_ = 0{N^), as pointed out in the 
calculation of the one-point Green's function. When we compute the following two level correlators: 

/— Tr i — Tr ^ \, (53) 

\N zi-H + ieN Z2-H-ie/' ^ ' 

the coefficient of Gaussian fluctuations become 1 — 2u^v^. The result of the saddle-point method 
should be trusted only for zi — Z2 = 0{N^) indicated by the deflnition of u± eq. (pO]) . Therefore, the 
short- distance two-level correlators should be calculated in other method instead of the saddle point 
evaluation. 



3 Gaussian unitary ensemble 

We turn to the case of the GUE, which means the ensemble of Hermitian matrices following the 
Gaussian distribution (|l|). The measure DH is explicitly written as 

DH = Y[d{ReHij)Y[d (Im Hij ) . (54) 

i<j i<j 
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We can proceed computation along the way in the case of GOE. We assume 



\z\ > 2 (55) 

in the case of the GUE, which corresponds to the assumption (|^). We can construct a generating 
function W{z) by employing complex vector variables as follows : 

W{z) ^/^1-Jd<P exp{-'-r^ {z - H) , (56) 

where 

n N 

nn^(Re#)rf(ImC)- (57) 

a=l i=l 

We can readily show that 

G{z) = hm{--^)^W{z) (58) 

in a similar way to the derivation of @. Difference between eq. (^ and eq. (^) arises because the 
complex vector variables contribute twice comparing to the case of real vector variables. We first 
take the average over H in (|56D and we introduce a Hermitian auxiliary matrix Qab along the line 
with the case of the GOE. See eq. ([10[) to eq. (|13|). Then the integration over the complex vector 
variables gives 

Wiz) = ^|DQexp|-^Tr(2 1og*(^-g) + g2)} 

= ^ Yld^aAiwf exp{--J2(2\ogt{z-Wa) + Wa')}. (59) 



Diagonalizing Q and we get the saddle-point equation corresponding to eq. (p!9|) 

1 



+ w; = 0. (60) 

z — w 

The most dominant saddle-point is w^In, in which 

U7_ = (61) 

We expand the exponent of ( ^9] ) around W-In and ignore the cubic and higher terms of the fluctuation 
as in the same way of GOE. The Gaussian integration of the fluctuation is easily performed, from 
which we derive 

W{z) ~ exp (21og(2; - w^) + wl)^ (j^^) (62) 



up to the next-to-leading order. Using this result and eq. (|58D , we get 



G{z) = ^{z-Vz^) + 0{1/N''). (63) 
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Note that first-order corrections vanish in this case. 

The two-point function defined in eq. ( pT]) is also calculated following the case of GOE. Using 
complex vectors, we define the generating function W{zi,Z2) as 

W{z„ z,) ^ l^Yiz; I ''''^^'2 '^^^ -H)<P1-\ 0f (^2 - H) <Pl}^ . (64) 

The two-point function G{zi, Z2) is derived from W{zi, z^) as 

We can write W{zi, z^) in terms of an integration over Hermitian (n + m) x (n + m) matrices, which 
corresponds to eq (EOl) for the GOE case. The result is 



The most dominant saddle point Q is 



W{z^, ^2) = ^ J DQ exp |-^Tr (2 \ogi{z - g) + g^) | . (qG) 



where 

= - - _ 4^ ^ ^ = 1, 2. (68) 

Since the connected part of the two-point function is the order 0(1/A^^), we must take in the effect of 
fluctuations around Q. As we explained in the GOE case, we can regard the fluctuations as Gaussian 
for the leading order of the connected part. Following the step to derive eq. (^Tp, we get 

/ \ \ ""^ 

W{z^,Z2) = W{z^)W{z2)[- . (69) 

VI - W1W2J 

Employing eq. (|65|) , the above result is translated to the language of the Green's function : 

G{zr,Z2) = G{z{)G{z2) - ^Q^^^og{l - G{z{)G{z2)) . (70) 
This agrees with the well-known result obtained by several other methods [Q, ^ ^ |rT|, |T^. 



4 Gaussian symplectic ensemble 
4.1 Definition 

The Gaussian symplectic ensemble(GSE) is the ensemble of quaternion real Hermitian matrices with 
the distribution (|l|). We shall first recall the definition of a quaternion real Hermitian matrix. 
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Let us write 



H 



( Hn 
H21 



HiN 

H2N 



(71) 



V Hni . . . Hnn / 

where each entry Hij is a quaternion number, which can be represented as a 2 x 2 matrix. We choose 
the set of the Pauh matrices and the unit matrix I2 as a basis of 2 x 2 matrices. Then Hij is 
written as 

3 



{k)k 



k=l 
rrll 7Tl2 

af3 



(72) 



In our notation, i, j = 1, - ■ ■ , N and a, /? = 1, 2 in H^f. The matrix H is called a quaternion real 



Hermitian matrix if and only if H^J^ forms a real symmetric matrix whereas H^j' {k = 1, 2, 3) form 
real antisymmetric matrices. These conditions are alternatively expressed as 



(k) 



H 



11* 



22 



Hi 



12* 



21 



Hji. 



(73) 



Note that ^ in the last equation means the Hermitian conjugation when we regard Hij cLS cL 2 X 2 
matrix. 

The trace of H'^ in eq.(|Tp reads 



(74) 



where a, (3 = 1, 2. Using -ffi^'', /i = 1, ■ ■ ■ 4, the measure DH in this case is written as 



i<j k=l i<j 



(75) 



Here, for later convenience, we consider how to construct a quaternion real Hermitian matrix 
from a complex vector with an even- number component. Let 0"" ^ i "£ N, 1 < a < n, a = 1,2 ) 
be a complex number. We regard 0"'s as the 2n-component vectors. The diadig 



or 



~a/3 



(76) 



defines the 2N x 2A^ Hermitian matrix. For constructing a quaternion real Hermitian matrix, we 
need to "symmetrize" k{(I)°') in the following way: 



n.22 

ij 



= K 



11 



11 



i fK21 _ S21\ ^ 21 



The matrix defined above is a quaternion real Hermitian matrix since it satisfies the condition 

^12* _ _ 21 

= Kji. (78) 

Another useful quaternion real Hermitian matrix is constructed as follows: we define 

iJi'''^<Pi'''*, (79) 

and make the diadig ipiipj. By the same symmetrization as eq. (|T^), we get the 2n x 2n quaternion 
real Hermitian matrix. We shall denote the resultant matrix by A(0i). A straightforward calculation 
gives 

Trft(0")2 = TrA(0,)2. (80) 

4.2 One-point function 

Let us compute the one-point function G{z) averaged over GSE 

G(z) = (^Tr^^). (81) 

Since if is a 2N x 2A^ matrix, Gi^z^ tends to I/2; as z ^ 00. We assume (|55|) as in the case of the 
GUE. The generating function W{z) is defined by 

W{z) ^/^1-Jd<P exp{-^0»t _ ^) (g2) 

Note that 0" is a 2A^-component vector, which implies 

(P''^ H (f)'' = (PI''* Htf <pf (83) 
in eq. (|S^) . The one-point function is derived from W{z) as 



Let us take the average over H in eq. (^2|) . We first notice that 

0'^t^0a = TrK(0")i7. (85) 
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We can symmetrize using the following relations 



which are derived from eqs. ([73D. For example, the coupling with the diagonal part of Hij in eq. (^) 
becomes 



'^"^ + ^I^f (87) 



After a similar procedure is performed to HlJ and Hfj, it turns out that the matrix /t(</)") defined 
in eq. ( [77| ) can replace in eq. (|85|). 

We are ready to take the average over H in eq .([821). That is, 

2 2 2 ^ 2A^ 8A^ ^ ^ 

The change H H + -^n{(f)°') can be performed by shifting the integration variables hI^'' because 
fi;(0'^) is a quaternion real Hermitian matrix. It leads to 

W{z) = Y^Jd<P exp [^-'-r^ z r - ^Tr K{<P'^f^ . (89) 

Next, for the sake of integration over 0, we introduce the auxiliary matrix Q that is a 2?t, x 2n 
quaternion real Hermitian matrix. Let us consider the following integration 



-^Jdq exp {-^q:m: + \i^rQ:^ih') , (90) 



where V'i's {i = 1, ■ ■ ■ n) are defined in eq. (^). 

The Q- integral in eq. (|90D is carried out in the same way as the if- integral of eq. (^21) . The result 

is 

^jDQ exp {-^Q:^QI: + \rrQ:^i^T) = exp (-^Tr \{rf) ■ (91) 



Thus, from eqs. (^,(|0l),(|9l) and (0), we conclude 

N 



W{z) = ^ J DQD<j> (-y Tr - '-^,\z - Q)^.) . (92) 
Performing the integration over 0, we get 

Wiz) = ^ //^Qexp|-yTr (2 1og^(z-Q) + Q2)}. (93) 
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Since any eigenvalue of Q has two-fold degeneracy 0, the eigenvalue representation of the above 
integration becomes 



1 „ n f n 1 

= ^ / W^Wa ^{wf exp {-NY, (2 \0gt{z - Wa) + Wa^) 
a=l I a=l J 



(94) 



Repeating the same argument as in the case of GUE, we obtain the result up to the next-to 
leading order: 

/ 1 \n^-n/2 

W{z) r~.e^^^^-Nn(2\og{z-w^)+wl)]{- — j , (95) 

where W- is defined in eq. (0). The exponent of (1 — w'^) is different from the result of GUE because 
the degree of freedom of fluctuations in the case of the GUE is given by ri^ while 2?t,^ — n in the case 
of the GSE. That affects the non- vanishing first-order correction to G{z). Namely, from eq. (p^), we 
obtain 



G{z) = -{z- .fz^A) (I - + 0{l/N^). (96) 

4.3 Two-point function 

Next we compute the two-point function G{zi,Z2) defined by 

The generating function W{zi, Z2) has the same form as the GUE case: 

W{z,, Z2) ^ / DhD<P2 exp{-'- 0f (^1 -H)^!-'- 0^ (^2 - H) . (98) 

The two-point function G{zi, Z2) is related to W{zi, Z2) as follows: 



In order to perform the integration over if, we introduce the notation as 



l<A<n 

lA—ti a 



^ bi-""" n + l<A<n + m ^ ' 



Then, since the coupling to H can be rewritten as we get the following result as the 

integration over H: 

W{z,, Z2) = / D<P,D<P2 exp{-^ z <l>^ - -^Tr ^i^^}, (101) 
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where z is the following 2(n + m) x 2(n + m) matrix 



Zil2m 



(102) 



According to the formulae ( ^Of ) and (0), we can express W{zi Z2) by 2(r2 + m) x 2(n + m) matrix Q: 



/% TV % 

D(I),D(I)2DQ exp{--<l>,t;2$, - -TrQ^ + -^^f^Q^,} 



1 r N i 

^— / D<P,D<P2DQ exp{--Trg2 _ ^^ _ g) v[/,}, 



where ^ is defined in the same way as eq. (ffOl): 



A2* 



After the integration over the vector variables, we get 

W{z,, ^2) = ^ / DQ exp |-^Tr (2 log2( 



Q) + Q' 



The most dominant saddle point Q is 



Q 



Wlhn 
W2l2m 



Next we compute effects of fluctuations around Q. 



W{ZI,Z2) = W{ZI)W{Z2) 



2nm 



, 1 — W1W2, 

Employing eq. (|65|), the above result is translated to the language of the Green's function 



G{z^,Z2) = G{z^)G{z2)- 



■\og{l-G{z,)G{z2)). 



(103) 



(104) 



(105) 



(106) 



(107) 



(108) 



2A^2 dzidz2 

This result with respect to the connected part agrees with that obtained by solving functional equa- 
tions 0, |ll| and by a diagrammatic method ||12[. Especially the disconnected part up to l/N"^ order 
in eq. (|108|) is consistent with [O]. 



5 Three-point function 

Finally, we calculate the connected three-point function in each Gaussian ensemble. That in the 
GOE case is defined by 

G(..,....3) - (>7^^T.^1t.^>. (109) 
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As is the case of one and two-point function, we construct the generating function of G{zi, Z2, Z3) 
using three n-fiavor real vectors: 

3 3 



W{zi,Z2,Z3) = 



/ n D<l>, exp{- ^ 1 ^^{zj - H)ct><^}). 

1=1 j=i 



The three-point function G{zi, Z2, Z3) is derived from W{zi, Z2, Z3) as 



_2 

G{zi,Z2,Z3) = hni(— ) 



^3 



-W{Zi,Z2,Z3). 



(110) 



;iii) 



n^o Nn dzidz2dz3 

We can write W{zi, Z2, Z3) in the term of an integration over 3n x 3n real symmetric matrix. The 
result is 



Wizi,Z2,Z3) = ^ [ L>gexp{-^Tr(logi(z-Q)+g2)}. 
The most dominant saddle point Q is 



(112) 



Q 



/ Wjn \ 

W2ln 
V Wsin J 



(113) 



where 



lizi-^Zi^-2), i = l,2,3. 



Next we consider fluctuations around the saddle point Q. 

W{zi,Z2,Z3) = ^ / D{SQ)e-^'-^Z3Si^ 



where 



So 

S2 



^iY(iogz(z-g) + g2), 

^TT{{5Qy-{Q6Qf), 

jj=l a,b=l 



1 



-iv(g5g)', 



z > 3. 



(114) 



(115) 



(116) 



(117) 



(118) 



Since the connected part of three-point function is in the order we must calculate the following 

Gaussian integral with 5*4 and terms: 

-So r °° 1 °° 



Wiz„Z2,zs) = J D{SQ)e-'^{l-j:Si + -{Y.Sir -■■■}, 

Q 1=3 ' 1=3 



e-Hi-m) + ^{{Sl)) + 0{l/N')}, 



(119) 
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where 



{{■■■)) = ^ J D{5Q){- ■ ■)e- 

Q 



S2 



(120) 



Calculating eq. (|119|) , we have 



n 



where 



jp/ X _ ^12 X23 X31 

[Zi, Z2, Zs) - —- —- — 

i — A12 i — A23 i — A31 

-^31 X12 1 ^ X12 X23 



(121) 



Xo:: 



23 



31 



1 — X31 1 — Xi2 1 — Xii 1 — X12 1 — X23 1 — X22 1 — X23 1 — X31 1 — X 



33 



-, (122) 



Xij = 2wiWj = ^G{zi)G{zj). 



Substituting the above formula into eg . ( p. 1 1| ) , we can derive the three-point function 

4 



; Z2, Z3) 



dzidz2dz3 



-F{Zi,Z2,Z3). 



(123) 



(124) 



The connected part in this result agrees with that obtained by other methods |T^, |TT], \L2 



Next we turn to the connected three-point function in the GUE case along the way in the case 
of GOE. We construct a generating function W{zi, Z2, Z3) by employing three complex vectors as 
follows: 

W{z,,Z2,z,) = I ^ f l[D<p,exp{-j2'-<pf{z, - H)<P'^}). (125) 

^(/.l ^02^03 i=l j = l ^ 



The three-point function is derived from W{zi, Z2, z^) in the following formula: 



^(21,2:2,2:3) = lim(— ^ 



93 



"^0 Nn dzidz2dz3 
W{zi, Z2, Z3) is written in the form of Q-integral 



W{Zi,Z2,Z3). 



W{zi, Z2, ^3) = ^ / DQ exp{-^Tr(2 logz(z - Q) + Q^)}, 



where Q means 3nx3n hermitian matrix. The most dominant saddle point Q is 

Q 



(126) 



(127) 



/ Wjn \ 

W2ln 
V W^In ) 



(128) 
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where 



i = 1,2,3. 



The dominant contributions to W{zi, Z2, Zs) becomes 



n 



W{zi, Z2, zs) = j^Fi^i^ ^2, Z3) 



(129) 



(130) 



where 



F{zi,Z2,Z3 

4 



X12 X23 X' 



31 



1 — X12 1 — X23 1 — X31 
-^31 X12 1 , X12 X23 



1 _|_ X23 X31 



1 — 1 — X12 1 — X^ 1 1 — X-i 9. 1 — 1 — Xr)/} 1 — X'},^ 1 — J^si 1 — X: 



-11 



-12 



23 ± — y\22 



^23 -L ~ ^31 



'-33 



-, (131) 



Xij = WiWj = G{zi)G{zj). 
Thus, the three-point function in GUE is 



Gc{zi, Z2, Z3) 



-F{zi,Z2,Z3). 



(132) 



(133) 



m dzidz2dz3 

The connected part in this result agrees with that obtained by other methods |[l|, [l^, |Tl|, [12|. 

To end with, we calculate the connected three-point function in GSE case. That in this case is 
defined by 



1 11 
G{zi,Z2,Z3) = {—Tt rr— Tr- 



1 1 1 V 

Tr -). 



'2N ZI-H2N Z2-H2N z^-H 
The three-point function derived from generating function is as follows: 

-1 93 



G{zi,Z2,Z3) = lini( 



n^o 2Nn dzidz2dz3 
And the generating function is calculated as 



W{Zi,Z2,Z3). 



W{zuZ2,Z3) = ^ J DQexp{-^Tr(21og2(z- Q) + Q2)}, 

Q 



(134) 



(135) 



(136) 



where Q means 6nx6n quaternion real hermitian matrix. The most dominant saddle point Q is 



Q 



( Wihn \ 

W2l2n 
V W^hn ) 



(137) 



where 



'>^i = 2^Zi - yzi^ - 4), 



i = 1,2,3. 



:i38) 
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Performing 1 /A^-expansion for W{zi, Z2, Z3), 



where 



W{zi, Z2, Z3) 



4n^ 



■F(ZI,Z2, Z3) 



F{zi,Z2,Z3) 



4X12 X23 X31 

1 — Xi2 1 — X23 1 — X31 



X31 X12 1 ^ X12 X23 1 X23 X31 1 



1 — X31 1 — X12 1 — Xii 1 — Xi2 1 — X23 1 — X22 1 — X23 1 — X31 1 — X 



33 



(139) 



(140) 



Xij = WiWj = G{zi)G{zj). 



The three-point function becomes 

; ^2, ^3 J 



1 



Q3 



2X4 dzidz2dz: 



-F{zi,Z2,Z3). 



(141) 



(142) 



The connected part in this result agrees with that obtained by other methods |11, 12 



6 Concluding remarks 

We have calculated the averaged one point Green's functions, the wide connected two-level and three- 
level correlators in Gaussian orthogonal, unitary and symplectic random matrix ensembles by the 
replica method. The one-point Green's functions have been calculated to the next-to-leading order 
in the 1/X expansion in GOE and GSE. Our results are consistent with those by other methods 
IQ, |TT], |T^. We have notified that there are some regions of the spectral parameter z where the 
employed saddle point evaluation cannot work well for averaged Green's functions. In those regions, 
the fluctuation of the auxiliary variable Q becomes large and the higher orders in the saddle point 
expansion give the same contribution with the order 0{N^). To calculate the short distance corre- 
lator we have to employ other method for the integration instead of the saddle point evaluation. 



19 



References 

[1] J. Ambj0rn, J. Jurkiewicz and Y.M. Makeenko, Phys.Lett.B251 (1990)517. 



[2] E. Brezin and A. Zee, Nucl.Phys. B402(1993)613. 

[3] M. L. Mehta, Random Matrices (Academic, New York, 1991). 

[4] C. W. J. Beenakker Nucl.Phys B422(1994)515; 

A. M. S. Macedo and J. T. Chalker, Physical Review B 49 (1994)4695. 

[5] J. Ambj0rn, C. F. Kristjansen and Yu. Makeenko, Mod.Phys.A7(1992)3187; 
E. Brezin and A. Zee, Nucl.Phys. B441 (1995)409; 
J. D'anna, E. Brezin and A. Zee, Nucl. Phys.B443(1995)443; 
E. Brezin and S. Hikami, Nucl PhysB479(1996)697 

[6] J. Ambj0rn and G. Akemann, J.Phys.A29(1996)L555 
G. Akemann, Nucl. Phys. B482(1996)403 

[7] S. F. Edwards and P. W. Anderson, J. Phys. F5(1975)965; 
A. J. MaKane and M. Stone, Ann.Phys. 131(1981)36. 

[8] J. J. M. Verbaarschot and M. R. Zirnbauer, J. Phys. A18 (1985) 1093. 

[9] T. Guhr, J. Math.Phys.32(1991)336. 
[10] J. J. M. Verbaarschot, H. A. Weidenmiiller and M. Zirnbauer, Ann. Phys. (NY) 153(1984) 367; 

[11] C. Itoi, to appear in Nucl. Phys. B, Universal wide correlators in non-Gaussian orthogonal, 



unitary and symplectic random matrix ensembles" |cond-mat / 96 1 1 2 T^ , NUP-A-96-15 



[12] C. Itoi and Y. Sakamoto, ^^Diagrammatic method for wide two level correlators in Gaussian 
orthogonal, symplectic random matrix ensembles" |cond-mat / 9 702 1 56 , NUP-A-97-5. 



20 



